Complex quantum network model of energy transfer in photosynthetic complexes 
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The quantum network model with real variables is usually used to describe the excitation energy 
transfer (EET) in the Fenna-Matthews-Olson(FMO) complexes. In this paper we add the quantum 
phase factors to the hopping terms and find that the quantum phase factors play an important role 
in the EET. The quantum phase factors allow us to consider the space structure of the pigments. It 
is found that phase coherence within the complexes would allow quantum interference to affect the 
dynamics of the EET. There exist some optimal phase regions where the transfer efficiency takes its 
maxima, which indicates that when the pigments are optimally spaced, the exciton can pass through 
the FMO with perfect efficiency. Moreover, the optimal phase regions almost do not change with 
the environments. In addition, we find that the phase factors are useful in the EET just in the case 
of multiple-pathway. Therefore, we demonstrate that, the quantum phases may bring the other two 
factors, the optimal space of the pigments and multiple-pathway, together to contribute the EET in 
photosynthetic complexes with perfect efficiency. 

PACS numbers: 87.15.H-; 89.75.Hc; 05.60.Gg 
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I. INTRODUCTION 

Photosynthesis provides chemical energy for almost all 
life on Earth. The initial step of photosynthesis involves 
absorption of light by the so-called light-harvesting 
antennae complexes, and tunneling of the resulting 
electronic excitation to the photosynthetic reaction 
center. Recent work has reported that quantum theory 
governs the exciton transfer in some light-harvesting 
complexes that harness the absorbed energy with almost 
100% efficiency The experimental evidence 

showing long-lived quantum coherences in this energy 
transport in several photosynthetic light harvesting com- 
plexes suggests that coherence may play an important 
role in the function of these systems. These observations 
have generated considerable interest in understanding 
the possibly functional role of quantum coherence effects 
in the remarkably efficient excitation energy transfer in 
photosynthetic complexes. 

The experimental achievements have motivated a num- 
ber of theoretical works (9l-[27j that consider the pho- 
tosynthetic complex as a quantum system, and try to 
analyze the basic mechanisms that explain the phenom- 
ena observed in the experiments. A full quantum dy- 
namic framework becomes necessary for studying coher- 
ent energy transfer. Typical quantum theories are, the 
quantum network model [9[-[l2l I14H19 . |25| . |27j , the hier- 
archic equation d, EH US Hl| j and the generalized Bloch- 
Redfield [22| equation. These theories have successfully 
predicted the long-lived quantum coherent phenomenon. 
However, all these models cannot describe the space dis- 
tribution of the pigments. As we know, the space dis- 
tribution of the pigments is very important for exciton 
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energy transfer and the experimental evidence also shows 
that the optimized space distribution of the pigments is 
one factor for prefect energy transfer in light-harvesting 
complexes @, E| . Therefore, it is necessary to set up a 
quantum model in which the space distribution of the 
pigments is considered. 

In this paper, we set up a quantum network model 
by adding the quantum phase factors to the two-body 
interactions to describe the exciton (the energy carrier) 
transfer in the FMO complex. The quantum phases are 
determined by the spatial structure of the pigments in 
photosynthetic complexes, such as the length of the pig- 
ment, the barriers and the distance between pigments. 
The quantum network in the absence of the phase factors 
has been used to study the EET in the photosynthetic 
complexes [9l-fl2l [l^ - [l9| . Some interesting results, such 
as the noises may enhance the EET, the EET in a quan- 
tum model may be larger than that of a classical model, 
are obtained. Compared with those studies, we find that 
the newly added quantum phase factors play a key role 
in the EET and there exist optimized phases at which 
the transfer efficiency is maximal. Furthermore, we find 
that the phase factors affect the EET just in the form 
of the phase difference in a closed loop. Although there 
may be many phase factors in the coupling terms of the 
system Hamiltonian, only N p — 1 are independent vari- 
ables where N p is the number of the pathways. It means 
that multiple-pathway in FMO complexes is a necessary 
condition for the enhancement of the EET by the quan- 
tum phases. It provides a strong evidence to support the 
statement that the multiple energy delivery pathway is 
also an acceptable contributing factor for perfect energy 
transfer Therefore, we demonstrate that, the quan- 
tum phases may bring the other two factors, the optimal 
space of the pigments and multiple-pathway, together to 
contribute the EET in photosynthetic complexes with 
perfect efficiency. 

The paper is organized as follows. In sec. II, we present 
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the complex quantum network we used to study the EET 
in photosynthetic complexes. A particularly simple and 
illustrative example with three sites of complex network 
is presented in sec. Ill, where some main conclusions, 
such as the phases may play an important role in the EET 
and the optimized phases are insensitive to the environ- 
ments, are demonstrated in this very simple example. In 
sec. IV, we present a symmetric complex quantum net- 
work with N p pathways to show that the multi-pathway 
is also a contributing factor for the perfect efficiency of 
the EET. In sec. V, we investigate the EET in FMO 
with our complex quantum network model. Finally, the 
conclusions are presented in sec. VI. 

II. THE QUANTUM NETWORK MODEL WITH 
THE PHASES 

The system we consider is a quantum network of N 
connected sites (nodes), schematically shown in Fig.l. 
Each site is modelled here as a spin-1/2 particle and it 
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FIG. 1: (Color online) The schematic representation of the 
quantum network which are a collection of N connected sites 
(nodes). Each site is modelled as a spin-1/2 particle (qubit). 
The particles are interacting with each other (solid lines) in 
the quantum network and may suffer dissipative losses as well 
as dephasing. An excitation is initialled at site 1. The arrow 
between site N and sink denotes an irreversible transfer of 
excitations from site N to the sink. 

may support an excitation which can be exchanged be- 
tween lattice sites by hopping. The initial (input) state 
is an excitation state which describes an excitation local- 
ized at site 1 (or several sites). We are interested to the 
transfer rate that the excitation transfers from the input 
state to the sink. The quantum evolution of the network 
of N sites is usually described by a Hamiltonian of the 
form 

N 

H = Y, wt°i + E v n tf°7 + (!) 

where and a~ are the raising and lowering operators 
for site j, crJ=\j)(Q\ and aj =|0)(j|, where |0) repre- 
sents the zero exciton state of the system and \ j) denotes 
the excitation being at the site j. The site energy and 



two-body coupling strength are given by the real num- 
bers Cj and Vji, respectively. The quantum network [28| 
described in Eq.([T]) has been used to study the EET in 
photosynthetic complexes in many literature EH - 

In this paper, we add a quantum phase factor e~ % ^ il 
with 4>ji a real number to the hopping term between sites 
j and I. The phase factor is determined by the detailed 
structure of the quantum network. As for the FMO, the 
phases are related to the length of the pigments as well as 
the intrinsic features of the barriers between the adjacent 
pigments. In this case, the Hamiltonian (1) is replaced 
by 

N 

Compared with the quantum network with real variables 
in Eq([T]). This model can be named as a complex quan- 
tum network model. We will show that the quantum 
phase factors paly the fundamental role in energy trans- 
fer of the photosynthetic complexes. 

As usual, the environment is modelled by the two dis- 
tinct types of the noise processes [Toj - TT^ ] . The first one is 
a dissipative process that transfers the excitation energy 
in site j to the environment with rate Tj , which leads to 
energy loss. The second one is a pure dephasing process 
with rate -fj which destroys the phase coherence of any 
superposition state in the system. The dissipative and 
the pure dephasing processes are described, respectively, 
by the Lindblad super-operators [Iol - [l2l |. 

N 

L dtss (p) = ]T ^[-{afa^p} + 2arpaf], (3) 

3=1 

N 

Ldeph(p) = ^1j[-{ (T i (T T>P} + 2(T j P a j ~ a j}> ( 4 ) 
J = l 

where {A, B} is an anticommutator. The absorption of 
the energy from the site k to the sink (numbered s) is 
modeled by a Lindblad operator 

L s {p) = T s [2(jfa-pCT+(7- - {(T+cr^cr+cr^p}], (5) 

where T s is the trapping rate. This term describes the 
irreversible decay of the excitations to the sink. So the 
full time evolution of the density matrix p of the system 
is described by the master equation 

^ = -~[H,p] + L diss (p) + L deph (p) + L s (p). (6) 

The efficiency of EET is measured by the popul ation 
Psink transferred to the sink from the site k \tfMY2\ , 

P sink = 2T s p kk {t)dt. (7) 
Jo 
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III. BI- PATHWAY QUANTUM NETWORK 

To study the role of the phase in quantum network, 
a particularly simple and illustrative example shown in 
Fig. 2 is to study quantum transport in a system of three 
sites. 

The exciton is transferred from site 1 to site 3 through 



where the phase difference 4> — 4>i2 + 4>23 ~ 4>i3i G = 
{DVB 2 + ACV 2 )[DT(T + T S )B 2 + C{2TT S + 7 r s + 
3DT)V 2 }, A = 3r + T s + 3 7 , B = 2T + T s + 2 7 , 
C = 3r + T s + 2-y, and D = T + T s . Although there 
are three phase factors </>i2, <^23, and (^13 in this three- 
site network, it is notable that only the phase difference 
4> = 4>\2 + 4>23 — 013 is independent. It demonstrates the 
fact that only the phase difference accumulated in the 
two pathways ( the pathway 1 — > 2 — > 3 and pathway 
1 — > 3 ) affects the interference at site 3. 
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FIG. 2: (Color online) The schematic representation of a bi- 
pathway quantum network. The exciton is transferred from 
site 1 to site 3 through two pathways, and finally, is trapped 
by the sink with the rate T s . 

From Eq. (8), we can easily find that P S i n k always 
increases with V, while it always decreases with T. How- 
ever, P S i n k is not monotonic functions of 7 and T s . We 
plot the dependence of P S ink on the different parameters 
in Fig. 3. 

Figure 3(a) shows the relation between P S mk and 7 
for different values of 4>. For <f) = 0, there exists an opti- 
mal value of 7 at which P s i n k is maximal, which indicates 
that the dephasing from the noise may even facilitate the 
EET. Note that the similar conclusion is extensively re- 
ported in the previous works [HI [HI [H| • When the phase 
4> is considered, the phase can change the efficiency re- 
markably at low dephasing (purely quantum mechanical) , 
while the efficiency is not sensitive to the phase at large 
dephasing (quantum coherent destroyed) . Therefore, the 
phase in EET plays a key role at low dephasing. 

Figure 3(b) shows the dependence of the transfer effi- 
ciency Psink on the trapping rate T s for different values of 



two pathways, and finally, is trapped by the sink with the 
rate T s - The sites I, 2, and 3 are susceptible simultane- 
ously to the dissipative and the pure dephasing processes. 
The dynamics of the system can be described by Eqs. (2- 
7). If we choose e\ = e 2 = £3 = e, V± 2 = V23 = Via = V, 
Ti = T2 = T3 = r, 71 = 72 = 73 = 7, an analytical 
expression of P S ink can be obtained 



(8) 

I 

<j>. When T s is very small, the system couples weakly to 
the sink, few exciton can reach the sink and the efficiency 
tends to zero. When T s is too large, the trapping rate 
F s mismatches the transport rate of the exciton in the 
quantum network, thus the efficiency also goes to zero. 
Therefore, there exists an optimal value of T s at which 
the efficiency takes its maximal value. 

Figure 3(c) and (d) show the efficiency P S ink as a func- 
tion of the phase <fi for different values of 7 and T, re- 
spectively. It is found that there are two optimal values 
(about 7r/2 and 37r/2) of <f> at which P S i n k takes its ex- 
ternum value, especially, it reaches a maximum value at 
37r/2. The minimal value of P S ink appears at <p ~ > 77 
and 2-7T. Obviously, the maximal values of P S mk are due 
to the constructive interference, while its minimal values 
are due to the destructive interference. Interestingly, the 
phases which correspond to the extremal values of P S ink 
almost do not change with 7 and T. 

Therefore, we can conclude that quantum phase in 
the two-body couplings plays a key role in the EET, es- 
pecially at low dephasing. Remarkably, the optimized 
phases in quantum network are almost independent of 
the environments (7 and T). 



IV. MULTIPLE-PATHWAY QUANTUM 
NETWORK 

Since the number of the pathways is an important 
quantity in quantum network, it is necessary to inves- 
tigate the role of the multiple pathways on the efficiency 
of EET. For simplicity, we consider a symmetric com- 
plex quantum network including N p pathways shown in 
Fig. 4. The exciton is transferred from site / to site F 
through multiple pathways, and finally, is trapped by the 
sink. N p describes the number of the pathways between 
the sites I and F. The dynamics of the system can also 
be described by Eqs. (2-7). From Eqs. (2-7), we can 
obtain the efficiency P S ink for different number N p of the 
pathways. 

Figure 5 shows the dependence of the transfer effi- 
ciency Psink on the phase <f> for different number N p of the 
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FIG. 3: (Color online) The dependence of the transfer effi- 
ciency Psink on the different parameters, (a) P S i n k vs 7 for 
different values of (j>. (b) Psink vs F s for different values of 
cj>. (c) P S i n k vs cj> for different values of 7. (d) P a ink vs </> for 
different values of F. Unless otherwise noted, the parameters 
are V = 1/5, F s = 1/5, T = 1/100, and 7 = 1/100. 



pathways. For a symmetric network, we choose ej = e, 
Vij = V Fj = V, Tj = T, 7j = 7, j = 1,2,3...JV P . There 
are N p — 1 independent phases because of iV p — 1 closed 
loops in the system. For simplicity, we only vary <f>n 
(= <f>) and the other phases are set to zero. For single 
pathway (N p = 1), the efficiency P S ink is always equal to 
0.695, this is due to the fact that no quantum interfer- 
ence can occur in single pathway. For double pathways 




FIG. 4: (Color online) The schematic representation of the 
symmetric quantum network including N p pathways. An ex- 
citation is initialled at site / and be transferred to site F 
through N pathways. The arrow between site F and sink de- 
notes an irreversible transfer of excitations from site F to the 
sink. 
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FIG. 5: (Color online) Transfer efficiency P 3 mk as a func- 
tion of the phase difference <f> for different number N p of the 
pathways. The other parameters are F = 1/100, 7 = 1/100, 
r s = 1/5, and V = 1/5. 



(N p — 2), the quantum interference at site F occurs and 
the phases take effect. Due to the destructive interfer- 
ence, there exists a minimal value of the transfer effi- 
ciency at <j) — 7T. The efficiency P S mk takes its maximal 
value at (f> = or 2tt, where the constructive interference 
occurs. The efficiency of multiple pathways at construc- 
tive interference has an enhancement compared with the 
single pathway. As the number of the pathways increases, 
the effects of destructive interference on the efficiency de- 
crease. We here have assumed that all other phases in 
N p > 3 pathways are zero. If we further optimize those 
phases, the enhancement by the multiple-pathway are 
clearer. It supports the conclusion in quantum scattering 
model[24j where the resonance transport is enhanced in 
multiple-pathway. Therefore, the multiple-pathway can 
reduce the destructive interference and facilitate EET in 
quantum network. 
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V. EXCITATION ENERGY TRANSFER IN FMO 
COMPLEX 

The architecture of antenna light-harvesting complexes 
varies widely among photosynthetic organisms. A well- 
studied example is the water-soluble FMO complex of 
green sulfur bacteria. The FMO complex essentially 
acts as a molecular wire, transferring excitation en- 
ergy from the chlorosomes, which are the main light- 
harvesting antennae of green sulfur bacteria, to the 



membrane-embedded reaction center. The FMO is a 
trimer made of three identical subunits, each containing 
seven pigments [13]. Because the inter-subunit coupling is 
vanishingly small, we only consider the dynamics of the 
EET within one subunit. The subunit containing seven 
pigments shown in Fig. 6(a) can be modelled as a net- 
work of seven sites with site dependent coupling and site 
energies. We use the experimental Hamiltonian of FMO 
given in [2{|, and the matrix of the Hamiltonian takes 
the form 
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with units of cm -1 and a total offset of 12230cm -1 to 
set the lowest site energy to zero for convenience (This 
overall shift in energy does not affect the dynamics of 
the system). In units with h = 1, we note that the rate 
lps -1 = 5.3 cm -1 . By neglecting the couplings weaker 
than 15 cmT 1 (only bold entries in the Hamiltonian are 
considered) in this model Hamiltonian, the transport in 
an individual monomer of FMO can be mapped to a 
quantum network shown in Fig. 6(b). 




FIG. 6: (Color online) (a) The spatial structure of one 
monomeric subunit of the FMO complex [f| Each 
monomer has seven pigments labelled by ® — (7). The initial 
state is taken to be a superposition state located at the pig- 
ments 1 and 6. and pigment 3 is in the vicinity of the reaction 
center (RC). (b) The simplified network for the monomeric 
subunit of the FMO. The thickness of two-head arrow in- 
dicates the coupling strengths and only couplings above 15 
cm -1 are shown. The exciton is transferred from sites 1 and 
6 to site 3 through the network, and finally, trapped by the 
reaction center with the rate T s . 



However, the Hamiltonian in Eq.® may not be suf- 
ficient to describe the EET in the FMO. We here fo- 
cus on the possible effects of the newly added phase 
factors in the coupling terms. From Eq. ([2|) and Fig. 
6(b) we can find that there are eight phases, 0i2, 023, 
034, 045, 047 , 056 , 067 , 016, but only two phase dif- 
ferences, 01 = 061 + 012 + 023 - 067 - 074 ~ 043 and 

02 = 067 + 074 — 065 — 054, are independent since there 
are just two independent closed loops. Therefore, with- 
out loss of the generality, we vary the phases 0i2 and 
067 and the other phases are set to zero in our numerical 
simulations. 

The initial state for our simulation is a superposition 
state localized at pigments 1 and 6 which are close to 
the chlorosome antenna (donor). It can be written as 
|*(0)) = a|l)+/3|6) with \a\ 2 + \/3\ 2 = 1. The pigment 3 
is the main excitation donor to the reaction center. The 
energy tra pping r ate from pigment 3 to the center in the 
literature [1IH12| ranges from 1 ps _1 to 4 ps _1 . In our 
calculations, we chose T s — 20/1.88 cm -1 corresponding 
to about 2ps _1 . The measured lifetime of excitons is of 
the order of Ins which determines a dissipative decay rate 
of 0.5/188cm _1 . Unless otherwise noted, we choose T = 
0.5/188cm _1 and 7 = 0.01T in this paper and assume 
that r and 7 are the same for each site. From Eqs. (2- 
7), we can numerically obtain the efficiency P S ink of the 
EET in FMO complex for different cases. 

Figure 7(a) and (b) shows the phase dependent effi- 
ciency for different dephasing and dissipative rates with 
067 = and a = f) = It is found that there exist 
two optimal values of 0i2 (e. g. 0i2 ~ 7r/2 or 3rr/2) 
at which the transfer efficiency P S mk takes its maximal 
value. When the dephasing rate 7 or dissipative rate T 
varies, the shape of the curve in Fig. 7 almost does not 
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FIG. 7: (Color online) Transfer efficiency P S ink as a function 
of the phase 0i2- (a)For different dephasing rates 7 at T — 
0.5/188 cm ~ . (b) For different dissipative rates F at 7 = 
0.5/18800 cm -1 . The other parameters are F s = 20/1.88 
cm" 1 , ^67 = and a = /3 = -^2. 



change, which indicates an important feature that the op- 
timized phases are not sensitive to the environment. The 
phases denote the distance between the pigments and the 
barriers between pigments, which are determined actu- 
ally by the spatial distribution of seven pigments. When 
the pigments are optimally spaced, the exciton can pass 
through the FMO with optimal efficiency Therefore, the 
phases from the two-body interactions play a key role in 
energy transfer of the FMO complex. 

Figure 8 shows the dependence of the transfer efficiency 
Psink on the phases </>i2 and 067 for different initial states. 
To study the significance of the phases, we can define the 
the difference AP between the maximal efficiency P™ 1 ^? 
and the minimal efficiency P™», AP = - P™«. 

We find that AP = 0.2467 for Fig. 8 (a) (a = /3 = 
AP = 0.6459 for Fig. 8 (b) (a = 0, (3 = 1) and 
AP = 0.1432 for Fig. 8 (c) (a = 1, fi = 0). Obviously, 
the phase can cause a significant change in the efficiency 
and the change just slightly depends on the initial states. 
The role of the phase coherence is to overcome local en- 
ergetic traps and aid efficient trapping exciton energy by 
the pigments facing the reaction center. In this cases 
we can still find that there exist some optimal phase re- 
gions where the transfer efficiency takes its maxima. The 
optimal phase regions only slightly vary with the initial 
states. Therefore, we demonstrate that the phase plays 
a significant role in the EET and the optimal phase can 
facilitate the energy transfer in FMO complex. 




FIG. 8: (Color online) Transfer efficiency P S ink as a function 
of the phases 0i2 and 067 for different initial states. (a)a = 
P = ^j-. (b)a = and /3 = 1. (c)a = 1 and /? = 0. The other 
parameters are F a = 20/1.88 cm" 1 , 7 = 0.5/18800 cm -1 , 
and T = 0.5/188 cm" 1 . 



The experimental evidences [H, Q show that beside the 
optimal space distribution of the pigments, the multiple 
energy delivery pathway is another acceptable contribut- 
ing factor for perfect energy transfer. Therefore, it is nec- 
essary to study the role of the multiple pathways in our 
model. In Fig. 9, the efficiency of one pathway is com- 
pared with that of the multiple pathways. It is found 
that the transfer efficiency P S ink gets a significant en- 
hancement when the number of the pathways increases. 
The efficiency is only 0.6425 for single pathway, while it 
can reach 0.8 for multiple pathways. Remarkably, the 
efficiency in the full FMO (case III) even at destructive 
interference is larger than that in single pathway (case I) . 
This result supports the conclusion in quantum scatter- 
ing model[24| where the resonance transport is enhanced 
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FIG. 9: (Color online) Transfer efficiency P S i n k as a function 
of the phase <f>56 f° r different cases at a — and /3 = 1. (I) 
One pathway: ©— >©— >d)— >■©, namely, V12 = V23 = Vie = 
V &7 = V47 = 0. (If) Two pathways: ®->®->-@->-(3) and 
®^©^(D^@, namely, V12 = V23 = Vie = 0. (Iff) Three 
pathways: the full quantum network of the FMO with cj>ie = 
O.lvr. Y a = 20/1.88 cm" 1 , 7 = 0.5/18800 cm" 1 , and Y = 
0.5/188 cm -1 . The other phases are set to zero. 



in multiple-pathway. Therefore, we can conclude that 
multiple pathways can also facilitate EET in FMO com- 
plex which agrees with the experimental statement 0, Q ■ 



VI. CONCLUDING REMARKS 

In this paper, we have investigated the efficiency of the 
EET in FMO complexes by adding the quantum phase 
factors to the quantum network model. The phase de- 
scribes the length of the pigments, the distance and the 
barriers between pigments and is then determined by the 



space distribution of the pigments. We found that the 
optimal distribution of the pigments can lead to the high 
efficiency of the EET. Moreover, the optimal phase is not 
sensitive to the environments. If the distribution of the 
pigments is optimized, the efficiency always takes it max- 
imal value, which is indeed significant for high transfer 
efficiency. As we know, the biological system governed 
by Darwinian selection has the optimal structure, which 
can ensure that the quantum coherence occurs in the op- 
timal spatial distribution. In addition, we also find that 
the multiple-pathway can facilitate EET in FMO com- 
plex. Therefore, we can conclude from the studies of the 
complex quantum network model that, the optimal space 
distribution of the pigments, the multitude of energy de- 
livery pathways and the quantum effects, are combined 
together to contribute to the perfect energy transport in 
the FMO complexes. 

However, in this paper we just add the phase factors 
phenomenally to the two-body couplings in the Hamil- 
tonian. How to determine them in a microscopic theory 
or from the experimental measurements is an important 
open question which deserves further study. 

Though we only have studied the transport process 
through the FMO protein, the methods and conclusions 
can be extended to other photosynthctic light-harvesting 
complexes. Furthermore, understanding the mechanism 
of efficient energy transfer in natural light-harvesting sys- 
tems can help developing low-cost and highly- efficient 
man-made solar energy apparatus, including photovoltaic 
devices and artificial photosynthesis. 
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